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Abstract:- Power system architecture today are being made more complex and they require more sophisticated 
new technologies or new devices such as Flexible AC Transmission Systems (FACTS). The limiting factor of the 
stability and the determination of the transfer limits depend on the load-voltage characteristic since load relief due 
to the load-voltage dependency results in larger transfer limits. This paper deals with the overview of a control 
strategies for the security assessment of an interconnected power system coordinated with different loads which 
can be governed by FACTS especially Interline Power Flow Controller (IPFC) devices when the system is 
approaching an extreme emergency state. IPFC is employed to enhance power system stability in addition to their 
main function of power flow control. In this method, the island is prevented from the total loss of supply using few 
FACTS devices. The optimization process is carried out using bacterial foraging optimization algorithm. The 
optimized result exhibits tremendous improvement in the system performance, when the proposed scheme is 
adopted in a IEEE 14 bus test system. 

Keywords: -Flexible AC Transmission System (FACTS), Interline Power Flow Controller (IPFC), 
Voltage Dependent Load (VDL) with Dynamic Tap Changer, Exponential Recovery Loads (ERL), Dynamic 
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I. INTRODUCTION 

Power system dynamic behaviour analysis requires adequate modelling of various utilities with the consideration 
of their parameter values that ensure those models with replicate reality. In many cases, the level of approximation is 
determined by the nature of the study. Load modelling provides an example. It is common for the aggregate behaviour of 
loads to be represented by a Voltage Dependent Loads with Dynamic Tap Changer and Exponential Recovery. The loads 
models will be a gross approximation, which provides information about the complex composition of the loads. This 
deficiency is particularly evident in distribution systems that supply a significant motor load, as the Voltage Dependent 
Loads with Dynamic Tap Changer and the Exponential Recovery Load model cannot capture the delayed voltage recovery. 
Design procedure arising from system dynamic behaviour can also be thought of in an optimization framework. However, 
any the optimization formulation must capture the processes driving dynamics. This class of problem be optimized 
efficiently using advanced computational procedures is referred as optimization algorithm. The classification of models is a 
decision that should be made based on the knowledge of the actual system composition and the phenomena that are being 
studied. 

In a power system network, there are often multiple loads connected to a single bus, normally the power of 
individual load is not measured or not available, but the total power transmitted through the bus is measured. In these cases, 
the loads can be considered as one composite load, which consists of static loads and dynamic or nonlinear loads. In recent 
years, many techniques have been proposed to model such loads . However, most of them are based on an assumed load 
equation and the parameters of the equation are estimated through curve fitting. Because of the complexity of modern loads 
(for example, power electronics loads), the assumed models may not capture power, frequency and voltage phenomena 
simultaneously and accurately. It is necessary to investigate new load modelling techniques to establish accurate load models 
for power system stability analysis. When a power system is subjected to large disturbances, such as [1] simultaneous loss of 
several generating units or major transmission lines, and the vulnerability analysis indicates that the system is approaching a 
catastrophic failure and more efficient and fast control actions need to be taken to limit the extent of the disturbance. 
Even with the FACTS controllers like Static Var Compensators, Unified Power Flow Controllers, Interline Power Flow 
Controllers trade-offs provides improved damping by optimization or tuning the controller even for small-signal conditions. 
As a consequence, the performance during the transient period immediately following a large disturbance may be degraded. 
FACTS devices' output limiters attempt to balance these competing effects. It can be shown that the tuning of these limiter 
values can be formulated as a bacterial foraging optimization problem. In contrast to design, analysis of system dynamics is 
more aligned with understanding extremes of system behaviour. The method has two stages. The first stage is to choose the 
weakest bus. The second stage includes a restoration process. Load restoration depends on the droop characteristic of the 
generators and maximum power transfer capability of the utilities. 

This paper discusses a vision for the solution for the power quality by improving the power system restoration 
through improved monitoring and control. The demand for larger power transfers over longer-distances, insufficient 
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investment in the transmission system, exacerbated by continued load growth are the trust areas which require continuous 
reforms to ensure a quality power supply to the consumers. The basic restoration assessment for Voltage Dependent Loads 
with Dynamic Tap Changer and Exponential Recovery Loads model in the power system network has been carried out and 
the various control corrective actions using IPFC is considered for the power system security enhancement study. 

II. POWER SYSTEM SECURITY ASSESSMENT 

Due to the nature of the disturbance and the set up of the power system network, there are two main elements in 
the Power System Security Assessment namely Static Security Assessment and Dynamic Security Assessment. Static 
Security Assessment is usually performed prior to dynamic security assessment. If the analysis evaluates only the expected 
post disturbance equilibrium condition (steady-state operating point), then this is referred as Static Security Assessment 
(SSA). Static security is related to an equilibrium point of the system where voltage and thermal limits are observed. It 
neglects the transient behavior and any other time dependent variations caused by changes in load conditions [2]. If the 
analysis evaluates the transient performance of the system as it progresses after the disturbance is referred as Dynamic 
Security Assessment (DSA). Dynamic Security Assessment is an evaluation of the ability of a power system to withstand a 
defined set of contingencies and to survive in the transition to an acceptable steady- state condition. DSA refers to the 
analysis required to determine whether the power system can meet the specified reliability and security criteria in both 
transient and steady- state time frames for all credible contingencies. 

In the operating environment, a secure system is one in which operating criteria are given prior important at pre- 
and post contingency conditions. This implies that analyses must be performed to assess all aspects of security, including the 
thermal loading of system elements, voltage and frequency variations (both steady state and transient), and all forms of 
stability. If the system is secured, these oscillations will decay and has to be damped out eventually. Otherwise, the 
oscillation of the frequency and voltage will grow to the extent of even shutting down the generator. If the system 
experiences a major disturbance, the oscillation will keep growing to a significant magnitude. The stability is then measured 
based on the trajectories of the disturbed systems motion related to the region of attraction of the final equilibrium state. For 
such situations, the accurate modeling of nonlinear system and regressive analysis of the nonlinear system [3] are to be 
adopted. Due to the tremendous growth in the interconnected Power system network with more complicated load models, 
increase the possibility of disturbance occurrences and the propagation of the disturbances. So the concept of the preventive 
(normal), emergency, and restorative operating states and their associated controls are to be adopted effectively. 

2.1. Off-Line DSA 

In the off-line DSA analysis, detailed time-domain stability analysis is performed for all credible contingencies 
and a variety of operating conditions. In most cases, this the off-line Dynamic stability analysis is used to determine limits of 
power transfers across the important system interfaces. Since the analysis is performed off-line there is no severe restriction 
on computation time and therefore detailed analysis can be done for a wide range of conditions and contingencies [4] . These 
studies include numerical integration of the models for a certain proposed power transfer condition and for a list of 
contingencies typically defined by a faulted location and specified fault-clearing time 

The simulation results are analyzed to find if voltage transients are acceptable and to verify whether the transient 
stability is maintained during the specified fault-clearing time. If the results for one level of power transfer are acceptable for 
all credible contingencies, the level of proposed power transfer is increased and the analysis is repeated. This process 
continues until the level of power transfer reaches a point where the system cannot survive for all of the credible 
contingencies. The maximum allowable transfer level is then fixed at the last acceptable level, or reduced by some small 
amount to provide a margin that would account for changes in conditions when the actual limit is in force [5]. 

III. MATHEMATICAL MODEL OF A POWER SYSTEM 

The dynamic behavior of multi-machine power system is described by the detailed modeling of all the elements of 
the power system [5]. 

A d 2 S dS _ _ 
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The general relation for dynamic behavior of a multi-machine power system used for both dynamic and transient 
security assessment be represented as 
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Where i = 1 ,2, N g . M h D t are the inertia and damping constant of the i th generation; P mi mechanical 

input to the i th generator; Ej is the EMF behind X' di of the i th generator; G^ B tj are the Conductance and Susceptance of the 
admittance matrix of the of the reduced system, X' di is the transient reactance of the i th generator; N g is the number of 
synchronous generator in the system. Modification of equ (3) for dynamic security assessment could result as 

m * s = p.-p.=p-p 
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3. 1. Modelling of Dynamic Loads 

A load can be mathematically represented with the relationship between power and voltage where the power is 
either active or reactive and the output from the model. The voltage (magnitude and/or frequency) is the input to the model. 
The load can be a static or dynamic load or a combination of both. Load models are used for analyzing power system 
security problems, such as steady state stability, transient stability, long term stability and voltage control. According to the 
power voltage equation, power system loads are divided into constant-impedance, constant power and constant current loads. 
A considerable amount of loads in power systems are induction motors and as several induction motors are connected to a 
busbar and it may be modelled as a single motor, parameters of which are obtained from the parameters of all motors and is 
referred as dynamic load modelling. [6]. 

The dynamic load model describes the time dependence as well as the voltage dependence of the load. The 
characteristic of a bus load depends on the load composition, which means that the aggregated load characteristics for the 
bus load must be found. These load parameters can be derived by field measurement based method which is based on direct 
measurement at a bus, during system disturbances or planned system disturbances, where voltage, frequency, active power 
and reactive power are measured. Then a method, such as the Least Square Method is used to derive parameters to the 
aggregated load model. 

Field Measurement Based Method is a simple Straightaway approach to derive a model. Unless the load 
composition is analysed in detail and unless buses having loads of fairly different compositions are measured there will not 
be performance analysis of the system and can even be extrapolated to different conditions. Spontaneous load variations are 
included in the load model, especially during long term measurements. Traditionally, lumped feeder loads are represented as 
static composite load models on the basis of constant impedance, constant current or constant power contributions. Models 
resulted in a classical dynamic model that was however still treated as linear in its voltage and frequency dependence of 
active and reactive powers [7] 

IV. MATHEMATICAL MODEL OF DYNAMIC LOAD 

The General form of Load Modelling equation can be written as, [8] 
x = a(x,v) 

P H =b p (x 9 v) 

Qh =b p (x,v) 
P H = b p O(oo), v) 

Where b p x(oo) solve x = 0, i.e. a(x, v) = for static load characteristic and, 

P H =b p (x(0),v) (9) 

Where P H is a Active power consumption model, x (0) is the value of the state when the initial change occurs. 
Dynamic load model with exponential recovery had also been proposed in [8], a differential equation that defines the 
behaviour of simple dynamic load model based on the response to a voltage step is given by [9]. 
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T p P H +P H =P s (V) + k p (V)V 

Where 

K p (V) = T p P t (v)v 

This equation can be written in first-order form as 
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Where T p is a Active load recovery time, P t (v) is a Transient part of active power consumption, X P is a state 
variable, P s (V) is the static load function . Using this form of differential equation, for active as well as reactive powers, the 
system voltage can be determined after the disturbance [9] . 

4.1.1. Voltage Dependent Loads with Dynamic Tap Changer 

Voltage Dependent Loads with Dynamic Tap Changer are nonlinear load model which represents the power 
relationship to voltage as an exponential equation. The transformer model consists of an ideal circuit with tap ratio n, hence 
the voltage on the secondary winding is vs = v/n. The voltage control is obtained by means of a quasi-integral anti-windup 
regulator The load powers P H and Q H are preceded as negative power as these powers are absorbed from the bus, as follows: 

-J jf = ^ol in) , 13 s 

~Qh - Qoi. In) q 4 x 
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and the differential equation is 

Where K d is the Anti-windup regulator deviation K { is the Anti-windup regulator gain. The reference voltage sign 
is negative due to the characteristic of the stable equilibrium point. If voltage dependent loads with embedded dynamic tap 
changer are initialized after the power flow analysis, the powers p and q are computed based on the constant PQ load 
powers p L0 and q L0 as (17andl8) and the state variable n and the voltage reference v re are initialized as follows 



(16) 



-m = V 

yref = -l + K d^ ^ 

Where v is the initial voltage of the load bus The parameters of this model are y p , y q , and the values of the active 
and reactive power, Po and Q , at the initial conditions. Common values for the exponents of the model for different load 
components y p and j q are (0, 1, 2). Equations can be directly included in the formulation of power flow analysis; where the 
y p and y q the active and reactive power exponent. The units of P and Q depend on the status parameter k. If k=l, the Voltage 
Dependent Loads with Dynamic Tap Changer is initialized after the power flow analysis, and P and Q are in percentage of 
the PQ load power connected at the Voltage Dependent Loads with Dynamic Tap Changer bus. 

P d = %00 ^ (17) 

if 

Qp = /ioo q l (18) 

Where K P and K Q are Active and Reactive power rating of the loads 

4.2.1. Exponential recovery load 

A nonlinear load model which represents the power relationship to voltage as an exponential equation as 



* P = -i + p*-p t 

Where p s and p t are the static and transient real power absorptions, which depends on the load voltage 

-« - * (V 
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where yt is a static active power exponent, ys is a dynamic active power exponent P s is a static reactive power 
exponent, p t is a dynamic reactive power exponent, of the load models, V is the actual voltage and V is the nominal voltage. 
Note that constant power, constant current and constant impedance are special cases of the exponential model. 
Similar equation holds for the reactive power 



■ : - - 
-Qb= ?+<?t 



(21) 
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(22) 



Where q s is the Static reactive load power as a function of bus voltage magnitude and q s is the Dynamic reactive 
load power as a function of bus voltage magnitude. The power flow solution and the PQ load data are used for determining 
the value of p , Q and V In particular p and Q are determined as below. A PQ load is required to initialize the 
Exponential recovery load bus. 

P D = P AoO ^ (23) 

Gd = AoO Q L (24) 

The parameters of the load can be defined based on the PQ load powers P L0 and Q L0 : 
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in this case initial voltage V is also not known, thus the following equation is used. 

-P H =gv 2 +I p v + P m (25) 
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2 
+ V ■ ^m (26) 



-Q H =bv 2 +Iv + Q ] 



The parameters are constants and indicate the nominal power is divided into constant power, constant current and constant 
impedance [10] 

4.2. Identification of Model Parameters. 

p h 4+K p (v-i)](i-p drop ) + p dyn .(Gv 2 -i) (27) 

Qh =k + K q (V-l)](l-Q dmp ) + Q dyn .(B-y 2 -1) (28) 

Considering the model given by equation (27), the nonlinear relationship between the measured signals, active power P, 
voltage V at the load bus, the estimated conductance G, and the parameters K P , P dwp and P dyn can be simplified by re- 
parameterization, and the model can be written as a linear regression equation 



p H =[x(D+x(2).(y-i)]+p dyn .(G.y 2 -i) 



(29) 



Where 

x(2) = x(\).K p 

W / \ / p ^ 30) 

6, = (x(l).x(2).P dyn ) (3i) 

The Least Squares method has then been used for the identification [10]. The objective is to obtain the best 
estimates for the parameter vector P , which minimizes the difference between the estimated active power and the simulated 
one (as a quadratic criterion). With the given equation (27) the same procedure is applied for the parameter identification for 
the reactive load using equation (28). The augmented objective functions to be minimized using a least square criterion 
which is given by equation (31). The final parameters are determined directly from the expressions given in (29) and (30) 
The least squares method is used to minimize the function (30) and to obtain the best estimates for the parameter vector 0p. 

*^{£'p ) 2L 4 \ simulated \ K ' & p ) *~ measure d\f K ' & p )) 

k=1 (32) 

The same procedure is repeated for the reactive powers also. The optimum solution represents the of the global 
minimum of the objective function, i.e. the best estimates for the model. However, the nonlinear model parameters can be 
estimated accurately by an iterative approach whose algorithm is as follows: 

* An initial estimate Xo for the parameters is selected; 

* The best fit is then determined by using the initial estimate Xo; the best estimates are compared with the initial estimates, 
and it is determined whether the fit improves or not. The direction and magnitude of the adjustment depend on the fitting 
algorithm [10]. 

V. OPERATION AND MATHEMATICAL MODEL OF INTERLINE POWER FLOW 

CONTROLLER UNIT 

The development of FACTS -devices [11-12] in ensuring high reliability as well as high efficiency with the modern 
power electronic components has elaborated the usage in various applications in the power system network. Voltage Source 
Converters (VSC) provides a free controllable voltage in magnitude and phase due to pulse width modulation of the IGBTs 
or IGCTs. High modulation frequencies allow to get low harmonics in the output signal and even to compensate disturbances 
coming from the network. The disadvantage is that with an increasing switching frequency, the losses are increasing as well. 
Therefore special designs of the converters are required to compensate this. Among various FACTS devices like SVC, SSSC, 
TCPS, UPFC it has been proven that even for multiline control IPFC contributes a lot [13]. So in this case study IPFC has 
been considered. 
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Fig.l .Power- injection model of FACTS device for power flow 

5.1. Steady state representation and load flow model of Interline Power Flow Controller (IPFC) 

When the power flows of two lines to be controlled, an Interline Power Flow Controller (IPFC) which consists of 
two series VSCs whose DC capacitors are coupled can be used. This allows active power to circulate between the VSCs. 
Fig. 2 (a, b) shows the principle configuration of an IPFC. With this configuration two lines can be controlled simultaneously 
to optimize the network utilization. In general, due to its complex setup, specific application cases need to be identified 
[14-15]. 
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Fig.2. (a) Schematic representation of an IPFC with two voltage source converters (b) Power injection model of IPFC 
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In this circuit, 

the complex Controllable series injected voltage is V sem = V s 

Series transformer impedance Z sem = R sem + jX sem 

the complex bus voltage at buses s m and r m be V rm = V sm ^0 rm and 

V sm = V sm <£_0 sm Z tm = R tm + jX tm and B tm represent the line series impedance and line charging susceptance, respectively 

[16-17], while m is the line number (m = 1,2....). From Fig (4), for one of the lines, the relations can be derived as follows 
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Fig.3. Equivalent circuit of IPFC with two voltage source converters 
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V tm and I sm can be expressed according to V rm and I rm as 
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From eqn. (33-36) I sm and I rm can be expressed in terms of V sm V rm and V sem I sm and I 
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Where V sm I sm V rm I rm and V tm are the complex bus voltages and currents at the corresponding buses s m and r 
respectively. As IPFC neither absorbs nor injects active power with respect to the AC system [18], the active power 
exchange between or among the converters via the dc link is zero, and if the resistances of series transformers are neglected, 
the equation can be written as 

(39) 
Where 

"dm — J^eai ( v sem i sm j. 

P dm is a active power exchange on the DC link, V sem is the series injected voltage and I sm complex bus current at bus sm . 
Thus, the power balance equations are as follows 

P + P -P -P = 

gm inj,m im tm (A(X\ 

Q +Q.. -Q -Q, = 

zzsgm z~sinj ,m -^im z^tm (AA\ 

Where P gm and Q gm are generations active and reactive powers, P Zm and Q /m are load active and reactive powers. V tjn and Q f m 
are conventional transmitted active and reactive powers at the bus m = " /" and "/'. 

1. Real power loss 

This objective consists to minimise the real power loss P L in the transmission lines that can be expressed as, 

^L-i&=itk (42) 

where Ns number of buses; Vk <oc k and Vh < oc h respectively voltages at bus k and h, Y^ and 9^ respectively modulus and 
argument of the k h -t h element of the nodal admittance matrix Y. 

2. Voltage deviation 

This objective is to minimize the deviation in voltage magnitude at load buses that can be expressed as, 

v = Y NL Iv - v ref l 

where NL number of load buses; ref Vi prespecified reference value of the voltage magnitude at the i-th load bus, ref Vi is 
usually set to be 1.0 p.u. 

VI. BACTERIAL FORAGING OPTIMIZATION ALGORITHM 

For over the last five decades, optimization algorithms like Genetic Algorithms (GAs), Evolutionary Programming 
(EP), Evolutionary Strategies (ES) which had drawn their inspiration from evolution and natural genetics have been 
dominating the realm of optimization algorithms. Recently natural swarm inspired algorithms like Particle Swarm 
Optimization (PSO), Ant Colony Optimization (ACO) have found their way into this domain and proved their effectiveness. 
Bacterial Foraging Optimization Algorithm (BFOA) [21] was proposed by Passion is inspired by the social foraging 
behavior of Escherichia coll Application of group foraging strategy of a swarm of E-.coli bacteria in multi-optimal function 
optimization is the key idea of the new algorithm [22]. Bacteria search for nutrients in a manner to maximize energy 
obtained per unit time. Individual bacterium also communicates with others by sending signals. A bacterium takes foraging 
decisions after considering two previous factors. 

The process, in which a bacterium moves by taking small steps while searching for nutrients, is called chemo taxis 
and key idea of BFOA is mimicking chemotactic movement of virtual bacteria in the problem search space. The control of 
these bacteria that dictates how foraging should proceed and can be subdivided into four sections namely Chemotaxis, 
Swarming, Reproduction, Elimination and Dispersal. These operations among the bacteria are used for searching the total 
solution space. 

6.1. Chemotactic Step 

This process is achieved through swimming and tumbling via Flagella. Depending upon the rotation of Flagella in 
each bacterium, it decides whether it should move in a predefined direction (swimming) or altogether in different directions 
(tumbling). In BFO algorithm, one moving unit length with random directions represents "tumbling," and one moving unit 
length with the same direction relative to the final step represents "swimming." The chemotactic step consists of one 
tumbling along with another tumbling, or one tumbling along with one swimming. This movement is can be described as 

o , (j+u,i)=e>(j.u)+c(f) 0L= 

MOM*) (44) 
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Where C (i) denotes step size ; A(i) Random vector ; A (i)Transpose of vector A(i). 

6.2. Swarming Step 

A group of E.coli cells arrange themselves in a traveling ring by moving up the nutrient gradient when placed 
amidst a semisolid matrix with a single nutrient chemo-effecter. The cells when stimulated by a high level of succinct, 
release an attractant aspirate, which helps them to aggregate into groups and thus move as concentric patterns of swarms 
with high bacterial density. The mathematical representation for E.coli swarming can be represented by 
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(45) 



Where 

J cc - Relative distance of each bacterium from the fittest bacterium 

S - Number of bacteria 

P - Number of Parameters to be optimized 

6 m - Position of the fittest bacteria 

d a ttract> w attract> h re peien b are the coefficients representing the swarming behavior of the bacteria. 

6.3. Reproduction 

The least healthy bacteria eventually die while each of the healthier bacteria asexually split into two bacteria, 
which are then placed in the same location. This keeps the swarm size constant. For bacterial, a reproduction step takes place 
after all chemotactic steps. 

N c +\ 
J 1 health = / J{i,j,k,l) 

j=1 (46) 

For keep a constant population size, bacteria with the highest /health values die. The remaining bacteria are allowed 
to split into two bacteria in the same place. 

6.4. Elimination and Dispersal 

In the evolutionary process, elimination and dispersal events can occur such that bacteria in a region are killed or a 
group is dispersed into a new part of the environment due to some influence. They have the effect of possibly destroying 
chemotactic progress, but they also have the effect of assisting in chemotaxis, since dispersal may place bacteria near good 
food sources. In BFOA, bacteria are eliminated with a probability of P ed . In order to keeping the number of bacteria in the 
population constant, if a bacterium is eliminated; simply disperse one to a random location on the optimization domain 
[23-24]. 



Problem Formulation 

Min J = k 2 Pi + k\ pi + k 
Subject to: 

P <P <p o <o<o V < V < V 

1 min — ± — ± max ; >^min — 5s^ — i^max ; v mm — ' — r max 

Where k , k b k 2 are cost coefficients 



(47) 



71 



Security Assessment of an Interconnected Power System Considering Voltage Controller 




Initialize variables, Set all loop-counters 
and bacterium index i equal to 



Increase elimination - dispersion loop 
counter 1 = 1+1 



Stop 




Increase Reproduction loop 
counter k = k + 1 




Perform 
Elimination - 
dispersal 
(Fori =1,2,.. 
S, with 

probability P ed , 
eliminates and 
disperse one to a 
random location) 



Increase Chemotactic loop counter j =j + 1 



Perform 

Reproduction (by 
killing the worse 
half of the 
population with 
higher cumulative 
health and splitting 
the better half into 
two) 
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Y 




Set bacterium index / = i + 1 




Yes 



Compute J(ij,k,l), adding the cell to cell 
attractant effect to nutrient concentration 
and set / last = / (ij-l,k,l) 



Tumble (let the i-th bacterium take a step 
of height C(i) along tumble vector A(i)) 



Compute J(ij+l,k,l), with the cell-to-cell 
attractant effect 



Set swim counter 
m = 




No 



m = m + 1 



Set 
m = N< 




Yes 



Swim (let the i-th bacterium take a step of 
height C(i) along the direction of the same 
tumble vector A (i)) 



Fig. 4 Flowchart of the Bacterial Foraging Algorithm 
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Read the system data 



From the bus admittance matrix Ybus 



Assume 5- for i =2, 3, 4 ,n 

And V? for i = 2, 3 m for PQ bus 



I 



Set iteration count K = 




Find Pi K and Qi K for i=2, 3, 4....n with 
IPFC device 



K 



K 



FindAPi fori = 2,3 n and AQi for i= 2,3, m 

Find APi K , AQi K for power flows in IPFC connected 

buses 



K 



K 



K A ^ K 



Find maxi APi maxi AQi and max APy , AQij 



Find slack bus power and all line powers and line 
flows 



Print result 




Fig.5. Flowchart for the Proposed Method 
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Form the conventional Jacobin matrix 



I 



Modify Jacobin matrix for 
incorporating the IPFC parameter 



i 



Solve for ( Al V I ) / ( I VI ), A5i 



i 



Update the bus voltage and the IPFC output voltages 




Set voltages at limit value 



VII. RESULTS 

The performance analysis of a IEEE 14-bus, 5-generator system coordinated with different types of Dynamic load 
models without and with IFPC Unit were studied and the optimum utilization requirement with the IPFC device for each 
load was determined using BFO technique. In this case of study the buses 5 and 14 are connected with VDL and ERL Loads. 
The system under the study is shown in Fig 6 and 7. 



-J^O 




®tj*i g_ 



© 



o 



Transmission Line #'s 
Bus #'s 



Fig.6. IEEE 14-bus test system one line diagram 
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CD 
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Figu.7. IEEE 14-bus test system with IPFC Units 
Table.l. Power flow solution for IEEE 14 Bus systems with VDL with Dynamic tap changer Load in bus4, 5 and bus 14. 



Bus 

No. 


Voltage 
Magnitude 


Voltage 
Angle 


Real 
Power 


Reactive 
Power 


1 


1.0300 


0.0000 


2.3045 


-0.4347 


2 


1.0000 


-5.962 


0.1830 


0.6624 


3 


0.9800 


-14.773 


-0.9420 


0.3103 


4 


0.9602 


-11.694 


-0.4780 


0.0390 


5 


0.9616 


-10.009 


-0.0760 


-0.0160 


6 


1.0000 


-16.570 


-0.1208 


0.1122 


7 


0.9766 


-15.192 


0.0000 


0.0000 


8 


1.0000 


-15.191 


0.0000 


0.1328 


9 


0.9607 


-17.071 


-0.3093 


-0.1740 


10 


0.9595 


-17.309 


-0.0900 


-0.0580 


11 


0.9757 


-17.075 


-0.0350 


-0.0180 


12 


0.9822 


-17.549 


-0.0610 


-0.0160 


13 


0.9753 


-17.597 


-0.1350 


-0.0580 


14 


0.9474 


-18.488 


-0.1490 


-0.0500 



Table. 2. Power flow solution for IEEE 14 Bus systems with Exponential Recovery Load in bus 4, 5 and bus 14. 



Bus No. 


Voltage 
Magnitude 


Voltage 
Angle 


Real 
Power 


Reactive 
Power 


1 


1.0300 


0.0000 


2.3502 


-0.4376 


2 


1.0000 


-6.071 


0.1830 


0.6812 


3 


0.9800 


-14.955 


-0.9420 


0.3166 


4 


0.9591 


-11.915 


-0.4780 


0.0390 


5 


0.9606 


-10.222 


-0.0760 


-0.0160 


6 


1.0000 


-17.105 


-0.1454 


0.1281 


7 


0.9747 


-15.586 


0.0100 


0.0100 


8 


1.0000 


-15.890 


0.0000 


0.1435 


9 


0.9585 


-17.493 


-0.3150 


-0.1860 


10 


0.9577 


-17.750 


-0.0900 


-0.0580 


11 


0.9748 


-17.562 


-0.0350 


-0.0180 


12 


0.9821 


-18.079 


-0.0610 


-0.0160 


13 


0.9750 


-18.166 


-0.1350 


-0.0580 


14 


0.9460 


-18.956 


-0.1490 


0.0500 
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Table.3. Power flow solution for IEEE 14 Bus systems with Exponential Recovery Loads in Bus4, 5 and Bus 14 with IPFC 

Unit 



Bus 

No. 


IPFC 


IPFC 


IPFC 


V 

p.u 


Angle 


Power at the 
Bus 


V 

p.u 


Angle 


Power at the 
Bus 


V 

p.u 


Angle 


Power at the 
Bus 


MW 


MVAR 


MW 


MVAR 


MW 


MVAR 


1 


1.0300 


0.00 


2.3287 


-0.1365 


1.0300 


-0.00 


2.2847 


-0.1361 


1.0300 


0.00 


2.3336 


-0.1467 


2 


1.0000 


-5.52 


0.2475 


-0.2012 


1.0000 


-5.44 


0.2473 


-0.2391 


1.0000 


-5.55 


0.2476 


-0.2364 


3 


1.0000 


14.38 


0.9420 


0.2950 


1.0000 


14.12 


0.9420 


0.2955 


1.0000 


14.29 


0.9420 


0.2953 


4 


1.0000 


11.88 


0.4780 


0.0446 


1.0000 


11.47 


0.4780 


0.0442 


1.0000 


11.68 


0.4780 


0.0443 


5 


0.9905 


10.06 


0.0586 


0.0236 


0.9919 


-9.76 


0.0588 


0.0238 


0.9924 


10.02 


0.0586 


0.0237 


6 


1.0200 


16.51 


0.1454 


-0.0303 


1.0300 


16.05 


0.1208 


0.0193 


1.0400 


16.82 


0.1454 


0.0291 


7 


1.0007 


15.35 


0.0100 


-0.0100 


1.0026 


14.66 


0.0000 


0.0000 


1.0020 


14.89 


0.0100 


-0.0100 


8 


1.0000 


15.34 


0.0000 


-0.0042 


1.0000 


14.65 


0.0000 


-0.0146 


1.0000 


16.57 


0.0000 


-0.0114 


9 


0.9923 


17.15 


0.3150 


-0.1861 


1.9950 


16.38 


0.3093 


-0.1740 


0.9949 


16.91 


0.3150 


-0.1860 


10 


0.9893 


17.35 


0.0900 


-0.0581 


0.9933 


16.63 


0.0900 


-0.0580 


0.9951 


16.98 


0.0900 


-0.0580 


11 


1.0008 


17.07 


0.0350 


-0.0180 


1.0078 


16.47 


0.0350 


-0.0180 


1.0137 


16.52 


0.0350 


-0.0180 


12 


1.0064 


17.50 


0.0610 


-0.0160 


1.0163 


17.08 


0.0610 


-0.0i60 


1.0291 


17.96 


0.0611 


-0.0160 


13 


1.0030 


17.68 


0.1350 


0.0580 


1.0121 


17.32 


0.1350 


-0.0580 


1.0267 


18.39 


0.1351 


-0.0580 


14 


1.0000 


19.08 


0.0499 


-0.0498 


1.0000 


17.84 


0.1499 


-0.0499 


1.0000 


17.58 


0.1490 


-0.0495 



Table.4. Power flow solution for IEEE 14 Bus systems with VDL with Dynamic Tap Changer Loads in Bus 4, 5 and Bus 14 

with IPFC Unit 



Bus 

No 


IPFC 


IPFC 


IPFC 


V 

p.u 


Angle 


Power at the 
Bus 


V 

p.u 


Angle 


Power at the 
Bus 


V 

p.u 


Angle 


Power at the 
Bus 


MW 


MVAR 


MW 


MVAR 


MW 


MVAR 


1 


1.0300 


0.00 


2.2837 


-0.1303 


1.0300 


-0.00 


2.2847 


-0.1361 


1.0300 


0.00 


2.2885 


-0.1404 


2 


1.0000 


-5.41 


0.2472 


-0.2102 


1.0000 


-5.44 


0.2473 


-0.2391 


1.0000 


-5.44 


0.2473 


-0.2441 


3 


1.0000 


14.21 


0.9420 


0.2951 


1.0000 


14.12 


0.9420 


0.2955 


1.0000 


14.12 


0.9420 


0.2955 


4 


1.0000 


11.65 


0.4780 


0.0444 


1.0000 


11.47 


0.4780 


0.0442 


1.0000 


11.45 


0.4780 


0.0442 
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5 


1.0220 


-9.84 


0.0588 


0.0237 


0.9919 


-9.76 


0.0588 


0.0238 


0.9928 


-9.81 


0.0588 


0.0238 


6 


1.0200 


16.01 


0.1208 


-0.0382 


1.0300 


16.05 


0.1208 


0.0193 


1.0400 


16.32 


0.1208 


0.0200 


7 


1.0019 


14.96 


0.0000 


0.0000 


1.0026 


14.65 


0.0000 


0.0000 


1.0032 


14.52 


0.0000 


0.0000 


8 


1.0000 


14.95 


0.0000 


-0.0109 


1.0000 


14.65 


0.0000 


-0.0146 


1.0000 


14.51 


0.0000 


-0.0181 


9 


0.9937 


16.73 


0.3093 


-0.1740 


0.9950 


16.38 


0.3093 


-0.1740 


0.9962 


16.16 


0.3093 


-0.1740 


10 


0.9904 


16.92 


0.0900 


-0.0580 


0.9933 


16.63 


0.0900 


-0.0580 


0.9966 


16.49 


0.0900 


-0.0580 


11 


1.0013 


16.60 


0.0350 


-0.0180 


1.0078 


16.47 


0.0350 


-0.0180 


1.0142 


16.52 


0.0350 


-0.0180 


12 


1.0064 


17.00 


0.0610 


-0.0160 


1.0163 


17.08 


0.0610 


-0.0i60 


1.0292 


17.47 


0.0610 


-0.0160 


13 


1.0030 


17.18 


0.1350 


-0.0580 


1.0121 


17.32 


0.1350 


-0.0580 


1.0268 


17.90 


0.1350 


-0.0580 


14 


1.0000 


18.61 


0.0490 


-0.0487 


1.0000 


17.84 


0.1495 


-0.0490 


1.0000 


17.11 


0.1490 


-0.0495 



VIII. CONCLUSION 

This paper presents the coordinated emergency control with the usage of FACTS device especially IPFC. It has 
been found that with the IPFC controller, the risk of load shedding is considerably reduced and can easily be adopted for 
emergency control. Moreover the result indicate that this comparison method successfully prevent the system from blackout 
and restore the system faster. 



APPENDIX 



Table. 5 .Generator Data [25] 



Generator 

Bus No. 


1 


2 


3 


4 


5 


MVA 


615 


60 


60 


25 


25 


xi (p.u.) 


0.2396 


0.00 


0.00 


0.134 


0.134 


r a (p.u.) 


0.00 


0.0031 


0.0031 


0.0014 


0.0041 


*j(p.u.) 


0.8979 


1.05 


1.05 


1.25 


1.25 


*d(p.u.) 


0.2995 


0.1850 


0.1850 


0.232 


0.232 


x"d(V-u.) 


0.23 


0.13 


0.13 


0.12 


0.12 


T'do 


7.4 


6.1 


6.1 


4.75 


4.75 


1 do 


0.03 


0.04 


0.04 


0.06 


0.06 


x a (p.u.) 


0.646 


0.98 


0.98 


1.22 


1.22 


x' (p.u.) 


0.646 


0.36 


0.36 


0.715 


0.715 


X" a (p.u.) 


0.4 


0.13 


0.13 


0.12 


0.12 


T 

± go 


0.00 


0.3 


0.3 


1.5 


1.5 


*■ go 


0.033 


0.099 


0.099 


0.21 


0.21 


H 


5.148 


6.54 


6.54 


5.06 


5.06 


D 


2 


2 


2 


2 


2 



Table.6. Bus Data [25] 



Bu 

s 

No 

• 


P 

Gene 
rated 

(p.u.) 


Q 

Gene 
rated 
(p.u.) 


P 

Load 
(p.u.) 


Q 

Load 
(p.u.) 


Bu 

s 

Ty 

pe 

* 


Q 

Generat 
ed 

max. 
(p.u.) 


Q 

Genera 
ted 

min.(p. 
u.) 


1. 


2.32 


0.169 


0.00 


0.00 


2 


10.0 


-10.0 


2. 


0.4 


0.424 


0.217 



0.127 



1 


0.5 


-0.4 


3. 


0.00 


0.234 


0.942 



0.190 



2 


0.4 


0.00 


4. 


0.00 


0.00 


0.478 


0.039 


3 


0.00 


0.00 
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5. 


0.00 


0.122 


0.076 



0.016 



3 


0.00 


0.00 


6. 


0.00 


0.00 


0.112 



0.075 



2 


0.24 


-0.06 


7. 


0.00 


0.174 


0.00 


0.00 


3 


0.00 


0.00 


8. 


0.00 


0.00 


0.00 


0.00 


2 


0.24 


-0.06 


9. 


0.00 


0.00 


0.295 



0.166 



3 


0.00 


0.00 


10. 


0.00 


0.00 


0.090 



0.058 



3 


0.00 


0.00 


11. 


0.00 


0.00 


0.035 



0.018 



3 


0.00 


0.00 


12. 


0.00 


0.00 


0.061 



0.016 



3 


0.00 


0.00 


13. 


0.00 


0.00 


0.135 



0.058 



3 


0.00 


0.00 


14. 


0.00 


0.00 


0.149 



0.050 



3 


0.00 


0.00 



* Bus Type: 1) Swing bus, 2) Generator bus (PV bus) and 3) Load bus (PQ bus) 

Table.7. Line Data [25] 



Fro 

m 

Bus 


To 
Bus 


Resistanc 
e p.u.) 


Reactan 
ce (p.u.) 


Line 

charging 

(p.u.) 


Tap 
ratio 


1 


2 


0.01938 


0.05917 


0.0528 




1 


5 


0.5403 


0.22304 


0.0492 




2 


3 


0.04699 


0.19797 


0.0438 




2 


4 


0.05811 


0.17632 


0.0374 




2 


5 


0.5695 


0.17388 


0.034 




3 


4 


0.6701 


0.17103 


0.0346 




4 


5 


0.01335 


0.4211 


0.0128 




4 


7 


0.00 


0.20912 


0.00 


0.978 


4 


9 


0.00 


0.55618 


0.00 


0.969 


5 


6 


0.00 


0.25202 


0.00 


0.932 


6 


11 


0.099498 


0.1989 


0.00 




6 


12 


0.12291 


0.25581 


0.00 




6 


13 


0.06615 


0.13027 


0.00 




7 


8 


0.00 


0.17615 


0.00 




7 


9 


0.00 


0.11001 


0.00 




9 


10 


0.3181 


0.08450 


0.00 




9 


14 


0.12711 


0.27038 


0.00 




10 


11 


0.08205 


0.19207 


0.00 




12 


13 


0.22092 


0.19988 


0.00 




13 


14 


0.17093 


0.34802 


0.00 





Table.8. Control Parameters of the Bacterial Foraging Algorithm 



S.No. 


Parameters 


Values 


1 


Number of bacteria ,S 


45 


2 


Swimming Length, N s 


4 


3 


Number of chemotactic steps, N c 


95 


4 


Number of reproduction steps, N re 


4 


5 


Number of elimination-disperse events, N ed 


2 


6 


Elimination and dispersal Probability, P ed 


0.25 


7 


w attract 


0.05 


8 


W re pelent 


12 


9 


^repelent 


0.02 


10 


^attract 


0.02 


11 


The run-length unit (i.e., the size of the step taken in 
each run or tumble), C(i) 


0.1 
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Table.9.IPFCData[15] 



Parameters 


Values (p.u./deg) 


Complex control series injected voltage 

V se l 


0.015 


Complex control series injected voltage 

v se2 


0.0216 


Phase angle 5 sel 


349.662 


Phase angle 5 se2 


185.240 


DC link capacitance voltage(V dc ) 


2 


Qini (rad) 


1.0878 
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